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АВ5ТКАСТ 


Іп this thesis several incompressible oscillatory flow and 
flutter problems were investigated. First, a previously 
developed unsteady panel code was modified so that systematic 
comparisons with Theodorsen’s classical theory could be 
accomplished. It was found that the panel code is in excellent 
agreement with the Theodorsen results. Second, the panel code 
was applied to the analysis of bending-torsion flutter. Again, 
general agreement with Theodorsen’s flutter predictions was 
obtained. In the experimental part of the thesis two flow 
visualization experiments were performed. First, the vortical 
flow patterns generated by an airfoil executing harmonic 
plunge oscillations were visualized. In the second experiment, 
the interference effects between a stationary airfoil anda 


small vane executing plunge oscillations were explored. 
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I. INTRODUCTION 


A. GENERAL 

In this thesis, several numerical methods were used to 
analyze the flow about an airfoil performing unsteady motion 
in an inviscid incompressible fluid. First, the unsteady 
motion of a Single airfoil was analyzed after modifying the 
U2DIIF code [ref.2]. The primary purpose was to verify the 
code against the proven theory of Theodorsen for analyzing the 
phenomenon of flutter. To accomplish this the U2DIIF code was 
modified to calculate aerodynamic values over a range of 
reduced frequencies and then apply these values to the 
flutter analysis. 

Next, the propulsive effects of a plunging airfoil were 
verified through experimental methods using a low speed 
plexiglas wind tunnel. 

Finally, an exploratory test was conducted in the 
department's smoke tunnel to study the interaction between a 


plunging airfoil and a stationary large airfoil. 


B. SCOPE 

Chapter II contains the modification of the single airfoil 
U2DIIF code into the code UPOT.f and extensive verification of 
this code against results produced by Theodorsen. Chapter III 


describes the UPOT code and explains the modifications which 


were added to solve the flutter determinant. In chapter IV the 
flow visualization experiment is described which was performed 
to study the vortical wake patterns produced by a plunging 
airfoil. In chapter V a second experiment is described which 
was performed to explore a plunging airfoil’s potential for 


control of flow separation. 


ІІ. SINGLE AIRFOIL ANALYSIS 


A. U2DIFF PANEL CODE 
1. Geometry 

Figure 2.1 shows a representation of the system that 
is analyzed using the panel code. Shown are the values for h 
(plunge) and o (AOA). 

2. U2DIFF 

The U2DIIF code was developed by TENG [ref.2] for the 
study of unsteady inviscid and incompressible flow over a 
Single airfoil. The code is based on the extension of the 
panel method, developed by Hess & Smith [ref.4] for steady 
potential flow problems, to include the unsteady motion of the 
airfoil that is continuously shedding vortices into the 
trailing wake. This vortex shedding process is nonlinear in 
that the wake vortices influence the flow over the airfoil 
which in turn alters the vortex shedding as the airfoil 
proceeds in time. 

The non-linearity of the unsteady flow makes this 
problem different from the steady flow problem which requires 
only simple Gaussian elimination. Teng developed a code that 
used an iterative type of solution. Typical program output 
includes the airfoil pressure distribution, force and moment 


coefficients, and the trailing vortex wake pattern. No 
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Figure 2.1 Airfoil Geometry 


attempt is made here to reproduce the work of Teng or to 
explore the operation of the U2DIIF code, but the reader is 


encouraged to review reference 2. 


B. PHASE PROGRAM 

The phase program was put through some verification by 
Neace [ref.9] and modified slightly in order to present 
results for harmonic motion. The code PHV3.f (phaseshift) was 
written by Neace to convert the time dependent output of lift 
and moment histories to harmonic output using an iterative 


curve fit algorithm: 


F(t)-Amp*Sin(ot-*4) (2.1) 


where Amp = amplitude of motion, o - frequency, and ¢ = phase 
angle between motion and the aerodynamic forces. One primary 
output of this program was the values of phaseshift (4) 
between the AOA and coefficients of lift (C ) and moment (C,) 
for the pitching airfoil and the phaseshift between the 
plunge value (h/2b) and the C and C, for the plunging airfoil. 
The other output was the amplitude O E the 


pitching or plunging case. 


C. MODIFICATION OF U2DIIF AND PHASE PROGRAM 
In an attempt to make the above mentioned codes more "user 
friendly", the two codes were combined into a single code 


named UPOT.f. The modification involved a new input file 


called UPOT.in which gives the user several options of 
operation. The input file can call for the analysis of steady 
flow only, straight and modified ramp motion, pitch 
oscillation, plunge oscillation, and the capability of 
performing the oscillation analysis over a series of reduced 
frequencies. A sample input file is shown in Figure (2.2). 

1. Output 

Outputs from the code have been limited to reduce the 

amount of computer space taken up by the code operation. A 
sample calculation was run using the input from Figure (2.2) 
and the on-screen output is shown in Figure (2.3). Of course, 
the user can modify the output portions of the code to 
minimize output. The following list describes the input/output 
files and the data they contain. 

a. UPOT.IN: The input file figure (2.2). 


D. CL.d: This file contains the various AOA values along 
with 1ts corresponding ee for each time Step. 


с. CM.D: This file contains the various AOA values along 
with its corresponding CG, for each time step. 


а. РНАЅЕ.а: This file contains the values of поп 
dimensional Came (с), AOA, -C Ci tor each time step. 

е. ЕОКО15.ОАТ: This file contains the values of non 
dimensional AOA, curve tut for Co quive fit for Сы, 
(used in the phase portion of program). 


f. CPSS.d: This file contains the steady state pressure 
coefficient for the mid point locations of all the air 
foil panels. 


а. CPU005.d: This file contains the unsteady pressure 
coefficient for the mid point location of all airfoil 


panels. (in this case the values are for an AOA equal 
to 5 degrees). 


PHZSWP.d: This file contains the phase information of 
the reduced frequency sweep portion for the program. 
Tic pr le contains pene, phase, angle өш с, апа су ава 
Бени Је НОВ И ico. 


FLUTTER.IN: This file contains information that can 
be used to solve the flutter determinant. It contains 
Kee Cee Cet ске, С Im TOF the pitch or plunge 
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ramp 
ramp 


RFREQ is basea on fuil chord 
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Sinusoidal pitch, 


n/a 
Translational harmonic oscillation 


Minimum/MAX AOA in degrees for IRAMP/ITRANS/IOSCIL 
MAX does not apply for ITRANS 


Location of Elastic Axis as a fraction of fuil chord 
Initial reduced frequency for program 


RFQSTP Reduced freq step size for a sweep of freq.'s(enter 0.0 if only one calcu 
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RFQFNL Final freq for the sweep 


DELHX 
DELHY 
DELI 

CYCLE ° 


NICICLE: 
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Translational] amount in the chordwise direction 
Max Translational amount in the vertical direction 
Min Translational amount іп the vertical direction (h/b) 


-In case of ramp, cycle=1.5 denotes airfoil 


(dist/full chord} 
(h/fullchord(b)) 


# of cycles for oscillatory motions 
is held 


at max aoa for the duration of .5 cycle 


-For steady state solution set it to 0 


$ of time steps for each cycie 
CYCLE*NTCYCLE is limited to 200 currently. 


Tolerance for convergence of the unsteady solution. (recommend using not 
Less Chan 
$ cf input aca for cp output 


.001) 


- angles should be in increasing order, 
- for oscilatory motions angles should increase 


first, then decrease. Decreasing angles are for 
the return cycle.. 





Figure 2.2 UPOT.IN 
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ATREOIL TYPE ° NACA 0012 АТВЕСТЕ 
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NO. PANELS UPPER SURFACE 
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FREQ SWEEP 
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Figure 2.3a 


0.993600 
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14 0183049593 12037547 -0.38682000584161 1.028793 
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.894883 
2413336 
4930159 
.945283 
.958651 
.970208 
„919910 
. 987718 
.993600 
2297525 
„999507 


0.893264 
0.894605 
9.896208 
0.898136 
0.900467 
0.903280 
0.906675 
0.910746 
0.915603 
9.921315] 
0.928123 
0.936086 
0.945516 
0.957020 
0.971841 


-0.086820 
-0.086820 
-0.086820 
-0.086820 
-0.086820 
-0.086820 
-0.086820 
-0.086820 
-0.086820 
-0.086820 
-0.086820 
-0.086820 
-0.086820 
-0.086820 
-0.086820 


05025328 
0.031258 
0.037745 
0.044881 
0.052175 
0.061567 
0.071411 
0.082565 
0.095332 
O OLS 
0.127830 
0.149558 
0.177698 
0.217400 
0.280753 


0.987255 
0.984247 
0.980946 
0.977302 
0973255 
0.968728 
0.963630 
0.957828 
0.951141 
0.943306 
0.933900 
0.922194 
0.906809 
0.884647 
0.848084 
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ово 


BESIN UNSTEADY FLOW SOLUTION 
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SLED 


— 
5 Оо р о о LD SA LU Ir“ 


— 
A 


alpha 


-3.0090 
=2 9836 
-2.9344 
“2.8532 
-2.7406 
22.9981 
=2 421] 
=2.2294 
-2.0074 
-1.7634 
-1.5000 
zI 2202 
=0.9271 
-0.6237 
-0.3136 
0.0000 
0.3136 
9. 6237 
359271 
1,2202 
1:5000 
1.7634 
2.0074 
2.2294 
224271 
2.5981 
2.7406 
2.8532 
2.9344 
2.9836 
3.0000 
2.9836 
2.9344 
2,8532 
2.7406 
259981 
2.8271 
2.2294 
2.0074 
1.7634 
1.5000 
122292 
029271 
0.6237 
0.3136 
0.0000 
-0.3136 
-0.6237 
-0.9270 
=]. 2202 
-1.5000 


time 


0.0000 
0.1518 
0.3035 
0.4553 
0.6071 
0.7588 
0.9106 
1.0624 
122141 
153659 
1.5117 
1.6694 
8212 
1.97130 
2.1241 
22165 
2.4283 
2.5801 
2.1318 
2.8836 
3,0354 
3.1871 
223289 
3.4907 
3.6424 
3.7942 
3.9460 
4.0977 
4.2495 
4.4013 
4.5530 
4.7048 
4.8566 
5.0083 
5.1601 
3. 3119 
5.4636 
5.6154 
5.7672 
5.9189 
6.0707 
6.2225 
6.3742 
6.5260 
6.6778 
6.8295 
6.9813 
91331 
7.2849 
7.4366 
7.5884 
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els 


-0.3479 
20:33 10 
-0.3243 
-0.3096 
2072929 
-0.2744 
=0-2542 
20:2323 
-0 2091 
-0.1848 
=0 1595 
-0.1334 
«0,1069 
-0.0802 
=0 20539 
-0.02 1 
-0.0012 
0.0238 
0.0478 
0.0704 
0.0916 
0.1110 
0.1285 
0.1439 
0:1571 
0.1679 
0.1763 
0,1822 
0,1855 
0.1862 
0.1843 
0.1799 
0.1729 
0.1636 
021520 
0.1383 
0.1226 
0.1050 
0.0859 
0.0654 
0.0437 
0,0211 
-0.0021 
-0.0257 
-0.0494 
-0.0729 
-0.0960 
“0. 1185 
=0. 1399 
-0.1602 
= 0:71 191 
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ста 


0.0002 
-0.0005 
20.0012 
-0.0020 
20:260217 
-0.0034 
- 22008» 
-0.0044 
-0.0047 
-0.0049 
-0.0049 
-0.0047 
-0.0043 
=0.0038 
=00031 
-0.0023 
-0.0014 
-0.0004 

0.0006 

0.0016 

0.0026 

0.0034 

0.0042 

0.0048 

0.0053 

0.0056 

0.0057 

0.0056 

0.0054 

0.0050 

0.0045 

0.0038 

0.0031 

0.0023 

0.0016 

0.0008 

0.0001 
-0.0005 
-0.0009 
=-0.0013 
-0.0014 
-0.0014 
=0.0013 
-0.0010 
-0.0006 

0.0000 

0.0007 

0.0014 

0.0021 

0.0029 

0.0036 


-0.0403 
=020425 
-0.0435 
-0.0442 
-0.0445 
-0.0446 
-0.0442 
2020436 
-0.0426 
-0.0412 
-0.0395 
-0.0375 
«020332 
-0.0326 
-0. 0291 
-0.0267 
-0.0234 
20.0199 
-0.0164 
-0:0121 
-0.0090 
-0:0052 
-0.0015 
0.0022 
0.0057 
0.0092 
0.0124 
0.0155 
0.0183 
0.0209 
9.0231 
020251 
0.0267 
0.0280 
0.0289 
0.0295 
0.0296 
0.0294 
0.0288 
9.0279 
0.0266 
0.0250 
0.0230 
0.0207 
0.0182 
0.0154 
0.0125 
020098 
0.0060 
0.0026 
-0.0009 





-1.7634 
-2.0074 
-2.2294 
-2.427C 
52.5981 
-2.7406 
228532 
-2.9344 


.7402 
.8919 
.0437 
22953 
„3472 
. 4990 
.6508 
.8025 


-0.0044 
0.0019 
-0.0113 
-0.014? 
-0.0179 
=00210 
= 0202 36 
-0.0265 


29.1963 0.0042 
=О0 2117 0.0047 
29.2251 (120051 
-0, 2 364 0.0053 
-0.2454 0.0054 
-0.2520 0.0053 
20.2562 0.0051 
-0,257В 0.0047 
0 


22.9836 
-3. 0000 
=> 38356 
-2.9344 
- 2. 8532 
-2.7406 
222288! 
“2,821. 
-2:2296 
z20014 
-1.7634 
-1.5000 
= 2202 
Ste EN 
-0.623! 
-0 3136 
0.0005 
0.3136 
06240 
0.9270 
12202 
1.5000 
1.7634 
2.0074 
2.2294 
2.4270 
25598; 
2.7406 
28552 
2.9344 
2.9836 
3.0000 
2.9836 
2.9344 
2.8532 
2.7406 
25981 
204211 
2.2294 
2.0074 
1.7634 
1.5000 
152202 
0.9271 
0.6237 
O. 3136 
0.0000 
EO. S16 
-0. 6237 
0. 9270 
-1,2202 
= 1 "5000 
-1.7634 
-2.0074 
-2.2294 
-2.4270 
-2.5981 
-2.7406 
22.823127 
-2.9344 
=29836 


„9543 
«1061 
„2578 
.4096 
5614 
a | 
.8649 
КОШЫ? 
. 1684 
42202 
.4120 
26233 
1/7 29 
49273 
11-0790 
11.8 
1722826 
2122242 
11.5851 
128479 
1.7891 
12.1414 
1272932 
12.4450 
12.9967 
12. 7485 
12.9903 
13,0220 
13.2038 
13.3556 
13.5073 
13.6291 
13.8102 
13.9626 
14.1144 
14.2662 
14.4179 
14.5697 
14.7215 
14.8732 
15.0250 
15:1768 
15. 3285 
15.4803 
15.6321 
15.7838 
195793926 
16.0874 
16.2391 
16.3909 
16.5427 
16.6945 
16.8462 
16.9980 
17.1498 
1723015 
17.4533 
17.6051 
17.7568 
17.9086 
18.0604 


bi + + — + 4 , — 
ОООООМЉО ФМ ч (O со со со (O оо 0 JA 
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оо 
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20.2570 
20.2536 
-0.2477 
-0.2394 
=0;2288 
зи 2158 
-0.2009 
-0.1839 
20.1652 
20.1439 
-0 1285 
-0, 1006 
-0.0770 
-0.0528 
-0.0282 
2020036 
0.0209 
0.0448 
0.0681 
0.0904 
051113 
051311 
0.1493 
0.1654 
uo. 1796 
DIST 
022012 
0.2085 
0.2133 
0.2156 
0.2153 
0.2125 
0.2072 
0.1995 
0.1894 
01770 
0.1625 
0.1461 
0.1219 
0.1081 
0.0870 
0.0647 
0.0415 
0:0] 78 
-0.0063 
2070305 
-0.0546 
-0.0782 
0. 1011 
720-1230 
-0.1431 
=0 21630 
=0: 180: 
-0.1965 
20.2103 
20 7209 
0.2312 
=052382 
-0,2421 
-0.2446 
-0.2441 
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% 
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Icd doas 
O O O O O O O O O O O O O O O O O о о о оо оо ос 


.0042 
. 0035 
.0028 
.0019 
- 0041 
· 0093 
‚0005 
„0012 
. 0018 
.0022 


0025 
0027 


.0026 


0024 


.0020 
- 0015 
. 0008 
. 0091 
. 0007 
.0016 
20025 
-0033 
- 0040 


0047 


-0052 
.0055 
2057 


0058 


.0056 
.0053 
.0048 
.0042 


0036 
0028 
0020 


.0012 
.0005 
.0002 
.0008 
.0013 
.0016 
.0018 


0018 


.0016 
. 0013 
.0008 


0002 


0.0005 
0.0012 
0.0020 
0.0028 
0.0035 
0.0042 
0.0048 
0.0052 
020055 
0.0056 
0.0055 
0.0053 
0.0049 
0.0044 


-0.0289 
20,0310 
2070328 
-0.0342 
-0.0354 
-0.0361 
2070365 
=) 0366 
-0.0362 
-0.0355 
20.03645 
-C.033C 
-0.0313 
20.0297 
2070269 
-0.0242 
-9. 0213 
ZU 5188 
"020557 
-0. 0116 
-0.0081 
-0.0045 
-0 0010 
0.0026 
0.0061 
0.0095 
0.0129 
0.0160 
0.0189 
00217 
0.0241 
0.0263 
0.0282 
0.0297 
0.0309 
0.0317 
020322 
0.0323 
0.0320 
0.0313 
0.0303 
0.0289 
0.0272 
0.0252 
0.0229 
0.0203 
0.0175 
0.0145 
0.0113 
0.0079 
0.0044 
0.0009 
-0.0026 
-0.0061 
2070096 
-0.0130 
*0.0163 
-0.0194 
7020223 
2070250 
-0.0274 





AMPLITUDE; clamp, cmarp : 0.2300964 3.4413978Е-02 
РНАЗЕ; cC lp; cmp з 91 - 371 85255 

AVERAGE DRAG, TOTAL DRAG : Lod Sor 2b=010 "9251697 966-02 
ETAS, WBAR : =0.2118947 77,2855391-09 
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2. Difference Between UPOT and U2DIIF/Phase 


The input file format was changed along with the 


following: 


e The program can now analyze a pitch, plunge or ramp motion 
that starts from any minimum value of Alpha or plunge 
(h/2b). Previously, the program only accepted the initial 
position of zero. This program does not need to go 
through the origin. 


e The phase portion of the program was changed to curve fit 
Сі апа С, Со а cosine function: 


F(t) =АтржСов (о Е+ф) (2.2) 


where, Amp = amplitude of motion, w = frequency of motion, 
ф = phase angle between motion and the aerodynamic 
forces. This was done since the alpha and plunge values 
were allowed to start from a new zero position. 


e The phase portion uses the middle 180 degrees of the final 
360 degree cycle specified in the UPOT.in file. This 
change was done to capture an all positive area of the 
cosine curve for phasing analysis. The program integrates 
this portion of the cosine curve, and for proper code 
operation the area under the curve must be kept to one 
sign. If the areas of integration were chosen to include 
both sides of the axis, then the code would produce errors 
near 90 and 270 degrees. 

3. UPOT Verification 
The code UPOT did not incorporate any drastic changes 
to the prior codes, but the original code had never been 
extensively compared to prior theories over a wide range of 
reduced frequencies. When conducting these comparisons, it is 


easy to become confused. This section will go through the 


comparisons slowly to help alleviate that problem. 
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а. Крапе1 (Кр) уз. Ктнеодогвоп (K+) 5 


The equation for reduced frequency is: 


_W2b _ wb 
аа Кр (2:31 


where: о = frequency of oscillation (rad/sec) 


b half chord (units to match U) 


U = free stream velocity (units to match b). 
The difference between Ko and Кү lies in the fact that Ko 
calls for the full chord and K, calls tor спе halt Ces 
hence, it is important to remember that Ko 1s twice K,. 

b. Aerodynamic Forces 

The aerodynamic forces problem of simple harmonic 

motion about an equilibrium position was solved theoretically 
by Theodorsen in NACA TR-496 [ref.10] and outlined by Fung in 
[геЁ.5]. The complex equations were simplified using the 


Simple harmonic motion equation and resulted in the following: 


L-npb?o? (Ls t [Lis (= +a) L,)a+[L,-(c-e) L,] B) ong (2.4) 
M-npb*«*( [M SM ] A. 
E i» 
[M,- (2 *a) (г, М,) + (3 «а) 21. | а (2.5) 


“[Mg- (S sa) -(с-е) ма (с-е) (S +a) L,) В) د ے‎ ) ۵٤ 


L, M are the lift and moment per unit span of the 


airfoil about the elastic axis, b, h/b, a and œ (radians), are 
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ое Dua Шш, La, and M, are defined by 
Scanlam зет 6 рр-а12 424) Бов various values of к. апа е 
This analysis will not cover airfoil aileron combinations. 


Therefore f becomes zero and equations 2.4 and 2.5 reduce to: 


L-npb?o? (L2 « (L,- (3 *a) L,] «) ei (e $0 (2.6) 
M=npb4w? ( [M,-(++a) L,] 2+ 
1 d В i (2.7) 
ІМ,-(--ға) (1,%М,) а) гр ја) е 1) 


The UPOT panel code used the following equations 


in defining lift and moment: 





T 
Ce -үР2 кі ейеефу Q stant (2.8) 
гаф Е. 
E M D 2 2 сі(шезфт) x -1 m 2 
C= – =у Ка +Тле ф„:Сап — (2.9) 
дар Ki 


here Road п are che real ang imaginary parts of C p and Ry 
and I, are the real and imaginary parts of C,. For the same 
conditions the lift (L) and moment (M) should be the same for 
both the panel code and Theodorsen. This fact allows for 
comparison of the magnitude, real, imaginary and phase of lift 
and moment. 


For litt: ре (eqn 2-6) equals Le (eqn 2.8) 


ES 


npb?w?[L,h/b+ (L,- (Veta) L,) a) e7 °°" saqb/Ri+Tpe* O°! (2.10) 


After canceling ei(9t*9D). 


npb?o? [L, (h/b) * (L,- (Ve*a) L,) «] -2qb| R2 * 12 


(2.11) 
бүк ет; 
For pitch case, h/b = 0, 2.11 reduces to: 
NPB O (AL - (Юта) Г.) а -2аБс, (2.12) 
Substitute К, = bw/U for e? and q = %pU? into 
equation 2.12 which gives: 
20 QbK; (L,~ (Veta) Г.) а-2аБС а (2.13) 
After cancelling and substituting Ky for Ee 
x K 
— P [L,- (Vata) L,] @=C,, (2.14) 


4 


This relationship can be further broken down into the real and 


imaginary parts: 





Ke 
Imag: T^P LiL, - (Veta) iL,] a=C..sin(o,) (2.15) 
4 а h La L 


Plunge Case: a=0 using equation 2.11 gives: 


20 


2 
x K. 
4 


Real: 





Іг,-(ға) І,|а-С,,со5(%,) (2.16) 


r p b362 L, =2qbC,, (2.17) 


The panel code uses h/2b for analysis because it uses full 


chord vice half chord. Therefore equation 2.17 becomes: 


2npb*o*L, (27) -2qbC,, (2.18) 


Substituting as before for o: 
2 fas 
2n qbKrL,2 (77) =2 GOC )2.19( 
Cancel and substitute Кр: 


Se? 


азе (2.20) 
| 2 ap! Do Сев 





This can also be broken up into imaginary and real parts as 
before. 


MOMENT : 


Equating equations 2.7 and 2.9 results in: 


npb*o*([M,- (Ve*a) L,] 23 ГИ а аео ва) 


=4 qb? ка + Та 


(2.21) 


For piten: h/D 0 
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8 EE S) 
м-4аЬ? ва + трае! (99 


си =R +12 


resulting in: 


праро ГИ ara (ЕМ) (бета) T Ju) it 9 


After substituting and cancelling: 


2 
атп Кр 


- Pe (Veta (eM) а у je 





«л К? 
8 





ІМАС: М, - 0 


2 
«т КЪ 





For plunge, @ = 0, equation 21 reduces го 


A 
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= [м„- (Уг+а) (1, +У2) + (а) 2Г.,| -С„соз (Фу) 


LIM ta ов, еее ЕС сш 


(2.22) 


(2.23) 


(2.24) 


(2.25) 


(2.26) 


(2.27) 


(2.28) 


IMAG: М, - 0 


ша A) (Weta) L,=CySin (yp) (2.29) 


2b 
Comparisons are shown for various cases of pitch and plunge. 
The tables include pitch values of 1 (Tables 2.2-2.3) and 6.7 
degrees (Tables 2.4-2.5), plunge (h/2b) values of .01 (Table 
2.8) and .0833 (Table 2.6-2.7). The graphs include 1 degree 
pitch (Figures 2.7-12, 2.13-20), 6.7 degree pitch (Figures 
2.21-28,2.29-36), .01 h/2b plunge (Figures 2.53-56), and .0833 
h/2b plunge (Figures 2.37-44, 2.45-52). 
4. Results 

The tables and graphs show that the panel code 
predicts the  Theodorsen results accurately. An initial 
question that was first addressed for the comparisons was how 
many cycles to use for good consistent phase results. Initial 
runs ere made at several different cycle values and the 
results are shown in Table 2.1. It can be seen that the panel 
code predicts Theodorsen’s results, using a cycle number of 
three. Increasing the cycle number just takes more computer 
time and only results in marginal increases in accuracy. 

The most glaring difference appears in the I, Cw 
comparisons of Figures 2.47-48. The panel code drops off 
sharply at the higher end of the reduced frequency spectrum. 
The reason for this is believed to be due to the magnitude of 


h/2b chosen for the comparison. The code was rerun for a 
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Comparison of Phase Calculations Using Various Cycles. 


| | їрисһ 6.7 Чер, МАСА 0007, .37c, 50 panels top and bottom! 
EEN | 


Ри Phase Angle |Cm Phase Angle |CÎ Amplitude |Cm amplitude 


54.409 04684 008927 

208.1592 -46.424| 0.5109224! 0.083932 
св: 206.001 0.51527| 0.083169] | |. 
ква 7 1.05% | 4.76% 0.84%) 0.92%, | 
| — i 


204.9365 -43.33887 0.51668, 0.082951 
204.3955 -42.868171 0.5172907 0.082865 
204.0596 -42.58497 0.517598! 0.082822; 


10. 0 203.8291 42.3916 0.517776) 0.082792 
as as 203.6031 -42.25684 0.51789} 0.082784 


% Diff. eon Te c. 0.11% 0.32% | 0.02% 0.01%. 


— 264.44052 -59.79002 1.0931! 0.212011 
3 261.9737 -58.29977 


0.94% | 22.56% 


260.17877 57.4404 1.10271; 0.21903 


259.2686 
258.5186 


257.9483 | 
257.6163 55.380567 1.106764) 0.218632 2 - | 


0.04% | 





TABLE 2.1 PHASE CALCULATION VS CYCLE NUMBER 


series of h/2b values and the percent difference for the panel 
code to Theodorsen was plotted in Figure 2.4. This chart shows 
that the h/2b value chosen has a tremendous impact on the code 
results. An h/2b value of .01 gave an acceptable error of 10% 
at K, = 8. Runs were completed with a value of .01 h/2b and 


the favorable results are shown in Figures 2.53-56. 
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зепел 32/4 





dy ЧбІҢ зе Аовапозу 
әрогу ио вәпірд GZ/y JO 12943 Әә) јо иозивдшо 


џезлоровц] 
шо 10ч84 10 НО 96 


Figure 2.4 Effect of h/2b values on Im C, 
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TABLE 2.2 1 DEGREE PITCH C, 
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Real Part of Cl for Panel and Theodorsen (pitch, 1.0 deg, .37c. 
МАСА0007, 75 panels top and bottom, 3cycles of 65 calc. per cycle) 
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Figure 2.5 1 Degree pitch С. Ке 
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Real Part of Ci for Panel and Theodorsen (pitch, 1.0 dag, .37c, 
МАСА0007, 75 рапаіз (ор and bottom, 3cycies of 65 calc. per сусів) 
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Figure 2.6 1 degree pitch C, Re 
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Imag. Part of Cl for Panel and Theordorsen (pitch, 1.0 deg, .37c, 


NACA0007, 75 panels top and bottom, 3cycles of 65 calc. per cycle) 
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Figure 2.7 1 Degree pitch C, Im 
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Imag. Part of CI for Panel and Theordorsen (pitch, 1.0 deg, .37c, 
МАСА0007, 75 panels top and bottom, 3cycles of 65 calc. per cycle) 
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Figure 2.8 1 Degree pitch C, Im 
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Mag. of Cl for Panel and Theordorsen (pitch, 1.0 deg, .37c, 
NACA0007, 75 panels top and bottom, 3cycles of 65 calc. per cycle) 
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Figure 2.9 1 Degree pitch C, Magnitude 
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Mag. of Cl for Panel and Theordorsen (pitch, 1.0 deg, .37c 
NACA0007, 75 panels top and bottom, 3cycles of 65 calc. per cycle) 
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Figure 2.10 1 Degree pitch C, Magnitude 
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Phase of Cl for Panel and Theordorsen (pitch, 1.0 deg, .37c, 
NACA0007, 75 panels top and bottom, 3cycles of 65 calc. per cycle) 
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Figure 2.11 1 Degree pitch C, Phase 
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Phase of CI for Panel and Theordorsen (pitch, 1.0 deg, .37c, 
NACA0007, 75 panels top and bottom, 3cycles of 65 calc. per cycle) 
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Figure 2.12 1 degree pitch C, Phase 
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TABLE 2.3 1 DEGREE PITCH C, COMPARISON 
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Cm Real vs Kp for Panel and Theodorsen (pitch, 1.0 deg, .37c, 
NACA0007, 75 panels top and bottom, 3cycles of 65 calc. per cycle) 
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Figure 2.13 1 Degree pitch Cy Re 
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Cm Real vs Kp for Panel and Theodorsen (pitch, 1.0 deg, .37c, 
NACA0007, 75 panels top and bottom, 3cycles of 65 calc. per cycle) 
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Figure 2.14 1 Degree pitch C, Re 
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Cm Imag vs Kp for Panel and Theordorsen (pitch, 1.0 
deg, .37c, NACAO007, 75 panels top and bottom, 
3cycies of 65 calc. per cycle) 
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Figure 2.15 1 Degree pitch C, Im 
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Cm Imag vs Kp for Panel and Theodorsen (pitch, 1.0 
deg, .37c, NACA0O007, 75 panels top and bottom, 
3cycles of 65 calc. per cycle) 
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Figure 2.16 1 Degree pitch С, Іш 
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Cm Magnitude vs Kp for Panel and Theordorsen 
(pitch, 1.0 deg, .37c, NACAQOO7, 75 panels top and 
bottom, 3cycles of 65 calc. per cycle) 
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Figure 2.17 1 Degree pitch C, Magnitude 
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Cm Magnitude vs Kp for Panel and Theodorsen (pitch, 
1.0 deg, .37c, NACAO007, 75 panels top and 
bottom, 3cycles of 65 calc. per cycle) 
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Figure 2.18 1 Degree pitch С, Magnitude 
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Cm Phase vs Kp for Panel and Theodorsen (pitch, 1.0 
deg, .37c, NACA0007, 75 panels top and bottom, 
3Jcycles of 65 caic. per cycle) 
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Figure 2.19 1 Degree pitch Cy Phase 
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Cm Phase vs Kp for Panel and Theodorsen (pitch, 1.0 
deg, .37c, МАСА0007, 75 panels top and bottom, 
3cycles of 65 calc. per cycle) 
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Figure 2.20 1 Degree pitch C, Phase 
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TABLE 2.4 6.7 DEGREE PITCH Ст COMPARISON 
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CI Real vs Кр Panel Code (Pitch, 6.7 deg, .37c, NACA 0007, 50 panels 
top and bottom) 
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Figure 2.21 6.7 degrees pitch C, Re 
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Cl Real vs Kp Panel Code (Pitch, 6.7 deg, .37, NACA 0007, 75 panels 
top and bottom) 
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Figure 2.22 6.7 degrees pitch C, Re 
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Cl Imag. vs Kp Panel Code (pitch, 6.7 deg, .37c, NACA 0007, 75 panel 
top and bottom) 
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Figure 2.23 6.7 degrees pitch C, Im 
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Cl Imag. vs Kp Panel Code (pitch, 6.7 deg, .37с, МАСА 0007, 75 panel 
top and bottom) 
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Figure 2.24 6.7 degrees pitch C, Im 
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CI Magnitude vs Kp Panel Code (pitch, 6.7 deg, .37c, NACA 0007, 50 
panels top and bottom) 
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Figure 2.25 6.7 degrees pitch C, Magnitude 
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Mag. of CL for Panel and Theodorsen (Pitch,6.7deg, .37 c, NACA 
0007, 75 panels top and bottom, 3cyc65calc.) 
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Figure 2.26 6.7 degrees pitch C, Magnitude 
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CI Phase vs Kp (pitch, 6.7 deg, .37c, NACA 0007, 75 panel top and 
bottom) 
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Figure 2.27 6.7 degrees pitch C, phase 
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СІ Phase vs Kp (pitch, 6.7 deg, .37c, NACA 0007, 50 panel top and 
bottom) 
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Figure 2.28 6.7 degrees pitch C, Phase 
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TABLE 2.5 6.7 DEGREES PITCH Cy 
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Cm Real vs Kp Panel Code (pitch, 6.7460, .37с, МАСА 0007, 50 
panels top and bottom) 


ТППЕ 
КШ 
КШ А 
ПИЯ 
ШШЩ 
ШЕШ 
ЕТШІ 


Kp = 2bw/U 


0.7 







0.6 


0.5 








0.4 





Cm Real 


0.3 
0.2 


0.1 


Figure 2.29 6.7 Degrees pitch C, Re 
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Cm Real vs Kp Panel Code (pitch, 6.7deg, .37¢, МАСА 0007, 50 
panels top and bottom) 
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Figure 2.30 6.7 Degrees pitch С, Ке 


55 


Cm Imaginary vs Kp Panel Code (pitch, 6.7deg, .37c, 
NACA 0007, 50 panels top and bottom) 
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Figure 2.31 6.7 Degrees pitch Cy Im 
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Cm Imaginary vs Kp Panel Code (pitch, 6.7deg,.37c, 
NACA 0007, 50 panels top and bottom) 
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Figure 2.32 6.7 Degrees pitch C, Im 
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Cm Magnitude vs Kp Panel Code (pitch, 6.7deg, .37c, 
NACA 0007, 50 panels top and bottom) 
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Figure 2.33 6.7 Degrees pitch Cy magnitude 
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Cm Magnitude vs Kp Panel Code (pitch, 6.7deg, .37c, 
NACA 0007, 50 panels top and bottom) 
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Figure 2.34 6.7 Degrees pitch C, magnitude 
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Cm Phase vs Kp Panel Code (pitch, 6.7deg, .37c, 
NACA 0007, 50 panels top and bottom) 
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Figure 2.35 6.7 Degrees pitch C, phase 
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Cm Phase vs Kp Panel Code (pitch, 6.7deg, .37c, 
NACA 0007, 50 panels top and bottom) 
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Figure 2.36 6.7 Degrees pitch C, phase 
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-.0833 C, COMPARISON 


6 PLUNGE h/2b 


TABLE 2 
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Real Part of CI for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.37 Plunge h/2b=.0833 C, Re 
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Real Part of Cl for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.38 Plunge h/2b=.0833 C, Re 
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Imag Part of Cl for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.39 Plunge h/2b=.0833 С, Іп 
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Imag Part of CI for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, МАСА0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.40 Plunge h/2b=.0833 C, Im 
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Magnitude of Ci for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.41 Plunge h/2b=.0833 C, Magnitude 
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Magnitude of Cl for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.42 Plunge h/2b=.0833 
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C, magnitude 
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Phase of CI for Panel and Theodorsen (plunge, .0833 
h/2b, .37c, NACA0007,75 panels top and bottom, 


3cycles of 65 calculations) 
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Figure 2.43 Plunge h/2b=.0833 C, phase 
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Phase of CI for Panel and Theodorsen (plunge, .0833 
h/2b, .37c, NACA0007,75 panels top and bottom, 
3cycles of 65 calculations) 
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Figure 2.44 Plunge h/2b=.0833 C; phase 
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TABLE 2.7 PLUNGE h/2b=.0833 См COMPARISON 


ZA 


Real Part of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.45 Plunge h/2b=.0833 C, Re 
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Real Part of СМ for Pane! and Theodorsen (plunge, 
.0833 h/2b, .37c, NACAO0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.46 Plunge h/2b=.0833 C, Re 
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тад. Part of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 








CM Imag. 








Figure 2.47 Plunge h/2b-.0833 C, Im 


74 








Imag. Part of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.48 Plunge h/2b=.0833 С, Іп 
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Magnitude of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.49 Plunge h/2b=.0833 C, Magnitude 
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Magnitude of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.50 Plunge h/2b=.0833 C, Magnitude 
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Phase of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.51 Plunge h/2b=.0833 C, Phase 


78 


Phase of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.52 Plunge h/2b=.0833 C, Phase 
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TABLE 2.8 PLUNGE h/2b=.01 См 
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Real Part of CM for Panel and Theodorsen (plunge, 
.01 h/2b, .37с, NACA0007,100 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.53 Plunge h/2b=.01 C, Re 
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Imag Part of CM for Panel and Theodorsen (plunge, .01 
h/2b, .37с, NACA0007,100 panels top and bottom, 
3cycles of 65 calculations) 
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Figure 2.54 Plunge h/2b=.01 C, Іш 
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Magnitude of CM for Panel and Theodorsen (plunge, .01 
h/2b, .37c, NACA0007,100 panels top and bottom, 
3cycles of 65 calculations) 

















0.08 
> 
© 
2 | 
т x 
o | 
0.06 | 
0.04 
А | 
0.02 
0 
0 1 2 3 4 5 6 7 
Kp = 2bw/U 


Figure 2.55 Plunge h/2b=.01 C, Magnitude 
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40 


CM Phase 
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Figure 


Phase of CM for Panel and Theodorsen (plunge, .01 
h/2b, .37c, NACA0007,100 panels top and bottom, 
3cycles of 65 calculations) 
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2.56 Plunge h/2b=.01 C, Phase 
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Comparison of CM Imag vs Kp for Various Values of Pitch 
(NACA0007,100 panels,1,3and 6.7 deg. pitch,3cyc65)(0c,) 
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Figure 2.57 Pitch 0c. 1.0,3.0,and 6.7 degrees vs Theodorsen 
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‘Comparison of Imaginary Dimensionless Aerodynamic Coefficient (Im 
My-8*Im Cmalpha/(alpha*pi*Kp^2)NACAO0007,100 panels,1,3,6.7 deg. 
pitch,3cyc65,0c.) 
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Figure 2.58 Dimensionless Aerodynamic Coefficient for 


1.0,3.0,and 6.7 degrees 
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ІІІ. FLUTTER DETERMINANT 


The proven accuracy of the UPOT Code enabled it to be used 


for the solution of the flutter determinant. 


А. FLUTTER THEORY 

In order to analyze the phenomenon of flutter, it is 
necessary to obtain the equations of motion of the system. To 
Simplify the problem the assumption is made that the actual 
motion of the system can be considered a combination of 
fundamental wing bending, and fundamental wing torsion. The 
system can then be replaced by an equivalent system containing 
an airfoil section of unit span restrained by springs against 
independent vertical motion (bending),and torsion as 
illustrated in Figure 3.1. This paper will not consider the 
aileron hinge case so f and c are set equal to zero. According 
to the class notes of M. Platzer [ref.1] the formulation 
proceeds as follows: 

Consider the balance of the elastic, inertial and 


aerodynamic forces on a mass element: 


e Total Inertial force: -/ dm(h"+ra")=- (Mh"+S.a") 
Mass:  M-Í dm 
Static Moment about the elastic axis: S. - | ram 


a 


e The moments about the elastic axis are: 
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-Ir(h"«ra")dm - -(r,e"*S h^) (3:58) 


Mass moment of inertia about elastic axis 


I, = {гат 


e Elastic restoring forces are:  -hC, zo 
The Equations of motion therefore become: 


I 


h"M+e"S,+hC, = Г 


| (322 
GIN SOC aí 


Where: C, = Torsional stiffness of the wing 
Cp = Stiffness of the wing in translation (plunge) 

M = Mass of the wing per unit span 
These equations can be written in a different way by 
expressing the spring constants in terms of the natural 
frequencies. Consider the airfoil to be so restrained that 


only one degree of freedom is permitted. The equations of 


motion become: 


C 
Mh'«hC, 0 во that Cor- E 
(3.3) 
// = h = С. 
1, Та С, < 0 зо thal о, = | — 
a 
Hence: Ch = Maj? C, = Тод“ 


The small structural damping of metal aircraft may be 
approximated by a force that opposes the motion and is in 


phase with the velocity. One assumes therefore that the 
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magnitude of the damping is proportional to the elastic 
restoring force. Since the motion of the airfoil is harmonic 
at the critical flutter condition, the structural damping can 
be accounted for by replacing the terms: 

ПЕ иЕе по (тар 

Ve wit GC (14 g.) 
Where g, and g, are damping constants multiplied by i to ensure 
that the damping force is in phase with the velocities in the 
simple harmonic motion. 
From equation 3.3 with: 

в) пе" тана (с) = беж 
Ме have: 

h" = оће!“ and а" « -огое!“" 


And the equations of motion become: 


eiet(-g?hM-o?aS,*hC,) =L (3.4) 


eiet(-g?aI -o?hS,*«C,) = M (3.5) 


The equations for the aerodynamic forces were given by 


Fung [ref.5] and are shown here: 


L-npb?o? (L2 « [L - (+a) Li] a) eiet (3.6) 


« 


Equating equation 3.4 to 3.6 and 3.5 to 3.7 yields: 
Substituting into equation 3.8 and 3.9 for C, and C, and 


using the following dimensional terms: 
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M-npb*o* ( [M,- ( = +a) Ly] => 
(327) 


[M, - (= +а) (L, *M,) “(жауы a) eitt 


а 


(-w2Mh-w20S,+hC,) =npb3w? (Ly > +[1,- (> +а) г] а) (3.8) 








(war ODS aC, *npb*o*^ ([M,- (= +a) EN ‚+ 
(3.9) 
Ім,-( ха) (L,*M,) + (> +а) г, а) 
р= М mE m. жа (3.10) 


i tt pb? “ Mb " N Mb? 


The equations simplify and after bringing all terms to the 


left result in: 


A—+Ba=0 
(3.11) 
р— +Еа =0 


This is a homogeneous equation whose solution is obtained 


if the flutter determinant is zero. 


А В 
- о (3.12) 
р Е 


Where: 


20 


2 
д-н(1-( 4) (94) (1512,1, 


& 


B=pX,+L,-L, (Vera) 
(3.13) 
D=px, +V¥2-L, (Vata) 


о) 


2 
E üz, [1=( SH (1*ig.)]-Ve(Vo*a) *M,-L, (Vae*ta) *L, (V2*a)* 





u is the ratio of the mass of the wing to the mass of a 
cylinder of air of a diameter equal to the chord of the wing. 


и and e, are the natural angular frequency (rad/sec) of 


a h 
torsional vibration around "a" (elastic axis) and the natural 
frequency in deflection, respectively. x, is the location of 
center of gravity of the wing measured from a. o is the 
circular frequency of wing vibration. 

The relationships between the code and Theodorsen derived 
earlier in Chapter II can be used here to simplify the 
equations: (note: no damping in this case g,-9g,-0) 


For A: manipulating equation 2.20, 


2С,, 


лк? ( D: (3.14) 


КРЕ 
resulting in: 


о) 2 2 2C 
даа OY c culter. 
а) о, лк? ( Ay 
P грб 





(3.15) 


For B: manipulating equation 2.14, 
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4c 
m (3.16) 
пКра 


L,-L,(V2ta) - 
resulting in: 


4C. 
дКра 


For D: manipulating equation 2.28 





4С 
-L, (1 =— MUR 
Le (asa) e Ay (3.18) 
т = 
о 
resulting in: 
- 4 Cun 
а en (3.19) 
nK (——) 
So 2p 
For E: manipulating equation 2.24 
2 8 Сма 
-Ya(VYat+a) +M,-L, (Veta) +L, (Vata) *=——> (3.20) 
amt Kp“ 
resulting in: 
O 8С 
Fore а = (3.21) 
0) алк? 


The determinant is expanded to AE-BD-0, and the real and 
imaginary parts are set equal to zero. Substituting (w/w) = 
X and solving the real (2 roots) and imaginary (1 root) 
equations for values of X corresponding to each reduced 


frequency value. These X values can be plotted as SQRT(X) 


- 


against K, and any intersections of real and imaginary parts 
Signify a flutter point. 


Knowing that 





K = 228 and ЖЕ =) (3.22) 


solve for seen 


2 Бо) 
О tical = —— (3.23) 
кух 


which is the critical flutter speed. 


В. UPOTFLUT CODE 
1. FORMULATION AND INPUT 

The equations derived in the flutter theory section 
above were programmed into a FORTRAN subroutine and attached 
to the UPOT.f code. The UPOT code was modified first to 
enable it to conduct a frequency sweep of pitch and plunge 
simultaneously. The resulting frequency sweep pitch and 
plunge array data is then sent to the flutter subroutine which 
provides the values of SQRT(x) and K, for plotting. The 


program also gives a best guess for the U based on the 


critical 
difference betweeen the real and imaginary SQRT (X) values. 
The input file UPOTFLUT.IN is very similar to the regular 
UPOT.IN file with the addition of actual physical properties 


of the system being analyzed. The user should start the 
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analysis in the pitch mode first (IOSCIL =1, ITRANS =0) to 
ensure complete coverage of all frequencies of interest. The 
following relations were taken from NACA TR-685 [ref.8] and 
should prove helpful in determining the physical properties 
needed for program operation. 

к - mass ratio - mpb?/M 


КЕ туу 


Outputs from the code have been limited to reduce the 
amount of computer space taken up by the code operation. A 
sample input and output file are contained in Figure 3.2. 
The following list describes the input/output files and the 


data they contain: 


a. UPOTFLUT.IN The input file figure 3.2a 


р. CL. ad Same as UPOT.f output 
се Mad Same as UPOT.f output 
d. PHASE.d Same as UPOT.f output 
e. CPSS.d Same as UPOT.f output 


tfo CBEOU5 d Same as UPOT.f output 
PHZSWP.d This file contains Ky, Pr, Pur Cra’ Cua 
PLHZSWP.d This file contains Кр, фт.Фы, Сьн, Смн 


i. PITCH.in This file contains Ko „Ста Ке, Сја Im, GO R E 
См Im 
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Aileron Torsion 
Spring at Aileron Hinge 


b = semichord (ft.) 
cb = distance between midchord and aileron hinge, positive if aft of 
midchord (ft.) 
eb = distance between midchord and aileron leading edge, positive 
Р - aft of midchord (ft.) 
ab = distance between rotation point (elastic axis) and midchord, 
positive if aft of midchord (ft.) | 
ћ = bending deflection of rotation point (elastic axis), positive down- 
ward (ft.) 
га = angular deflection about rotation point (elastic axis), positive 
for leading edge up (radians) 
В = angular deflection of aileron about aileron hinge relative tu wing 
chord, positive for aileron leading edge up (radians) 





Figure 3.1 Simplified System Geometry 
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оо ое р. 
AIRFOIL TYPE : NACA 0012 AIRFOIL 
NLOWER = 50 , NUPPER = 50 
$ à à à à 9 9 6 9 $ 9 $ à 9 à à à à à à € à à à à à à & à à à $ 0 é 00 4 é £ ç 0 9239.5 0 @e 0 e o ea 00 é é ç ç 0e@  é# à 0 0 0 é é# à ó 
IFLAG NLOWER NUPPER 
0 50 50 
AIRFOIL TYPE 
7 
IRAMP IOSCIL ALPI 
0 1 -3.0 
FREO RFQSTP RFOFNL 
‚85 201 „95 
1605Т UGUST  VGUST 
0 0. Ds 
ITRANS PDEL!X  DELHY DELI PHASE 
0 0.00 .0833 -.0833 0.00 
CYCLE MICYCLE TOL. 
3 65 0.005 
naot & naot X aoa values multiplied by 10 (integer) 
2 05. 10 20 25 39:50 
Semi-chord Walpha Wh Mass 
6 90 72.0 .53789 
Ialpha Salpha Density 
4.84102 - 645468 .002378 
Comments... 


I RAMP : п/а RFREQ is based on full chord 
Straight ramp 
Modified ramp 


IOSCIL n/a RFREQ is based on full chord 
Sinusoidal pitch, motion starts at min Aoa 


ITRANS : n/a 
Translational harmonic oscillation 


§ of cycles for oscillatory motions 
-In case of ramp, cycle=1.5 denotes airfoil is held 
at max aoa for the duration of .5 cycle 
-For steady state solution set it to 0 


NTCYCLE: f$ of time steps for each cycle 
CYCLE*NTCYCLE 1s limited to 200 currently. 


NAOT: | of input aoa for cp output 
- angles should be in increasing order, 
- for oscilatory motions angles should increase 
first, then decrease. Decreasing angles are for 
the return cycie.. 


SEMI-CHORD Half Chord in feet. 


Walpha,Wh, uncoupled natural frequencies of the system in question. 
Walpha is pitch and Wh із різпде (Н2). 


Mass specific mass of the system in slugs/foot of span 


Та! рна Moment of Inertia of system about the elastic axis (a) 
per unit span length. 
Sal pha Static moment of wing-aileron per unit span length 


- 


Density Mass of air per unit of volume(slugs per ft^3) 





Figure 3.2a UPOTFLUT.in example input file 
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ФеедедееееоееФ%Ф%е%%%%89е99ее888е999468%98999%9еее8ееедеегеееее%9%9е64е9292246 


NACA 0012 AIRFOIL 


AIRFOIL TYPE 
NLOWER 50 


= 
= 


, 


NUPPER = 50 


%дезе%%8%4ее999%4е%е8%егееееедееегееее%%Ф%ееде9%990%69е998098968е0009992% 


OSCILLATORY MOTION, 


FREQ SWEEP 
FREQ 0.850000 


IOSCIL 


1 


PHASE SHIFT ANALYSIS 


FREQ 
0.8500000 
0.8500000 


ы 
kp- 1 freq 
AMPLITUDE; 
foscil Oitra 
PHASE; clp, cmp 
AVERAGE DRAG, TOTAL DRAG 
ETAS, WBAR 


clamp, cmamp 


- 


ns 


0.8500000 
1 


0.2316691 


2.9221054E-02 


202.4033 
1.5930884Ғ-03 
-0.1430808 


-63.70797 
0.1051438 
-1.1134184Е-02 


PHASE SHIFT ANALYSIS 


FREQ 
0.8500000 
0.8500000 


ы 
Кр- {freq 
AMPLITUDE; 
loscil Oitra 
PHASE; cip, cmp 
AVERAGE DRAG, TOTAL DRAG 
ETAS, WBAR 


clamp, cmamp 


FREQ SWEEP 
FREQ 0.860000 


n 


3 


0.8500000 
1 


0.2779089 2.2117507Е-02 


1 
270.5030 37.10253 
-6.2051453Е-03 -0.4095396 
0.1050769 -0.1181067 


PHASE SHIFT ANALYSIS 


FREQ = 
0.8600000 
0.8600000 


ы 
Ер- 1 Ғтеа 
AMPLITUDE; 
ioscil Oitra 
PHASE; clp, стр 
AVERAGE DRAG, TOTAL DRAG 
ETAS, WBAR 


- 


clamp, cmamp 


- 


n 


s 


0.8600000 
2 
0.231267 2 9490597Е-02 
i 202.8760 х 


299216965203 
-0.1398697 


-63.85251 
0.1050829 
-1.1383206E-02 


PHASE SHIFT ANALYSIS 


FREQ = 
0.8600000 
0.8600000 


w 
kp= ifreq 
AMPLITUDE; clamp, cmamp 
ioscil Oitra 
PHASE; сір, стр 
AVERAGE DRAG, TOTAL DRAG 
ETAS, WBAR 


FREQ SWEEP 
FREQ 0.870000 


ns 


. 
° 
° 
° 
° 


0.8600000 
2 


0.2804866 2.2494521Е-02 


= 1 
270.7843 36.19393 
-6.3188463E-03 -0.4170439 


0.1047920 -0.1205979 


PHASE SHIFT ANALYSIS 


FRFQ 
0.8100000 


w 


0.8700000 





Figure 3.2b UPOTFLUT example output file 
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kp= 0.8700000 ifreq 
AMPLITUDE; clamp, cmamp 3; 
foscil = Oitrans 
PHASE; clp, стр : 
AVERAGE DRAG, TOTAL DRAG : 


ETAS, WBAR 


0.2317936 
= 0 
203.3604 
1.5893428Е-03 
-0.1365460 


PHASE SHIFT ANALYSIS 


FREQ = 
0.8700000 
0.8700000 


ы 

kp= ifreq 
AMPLITUDE; clamp, cmamp : 
foscil = Oitrans 
PHASE; clp, cmp 5 
AVERAGE DRAG, TOTAL DRAG > 
ETAS, WBAR 2 


FREO SWEEP 
FREQ = 0.880000 


0.8700000 


3 


0.2830637 


271.0655 
-6.4334869E-03 
0.1045149 


PHASE SHIFT ANALYSIS 


FREQ = 
w 
kpe 


0.8800000 

0.8800000 ifreq 
AMPLITUDE; 
foscil = Oltran 
PHASE; сір) стр 
AVERAGE DRAG, TOTAL DRAG 
ETAS, WBAR 


clamp, cmamp : 


0.8800000 
4 
0.2318744 


203.8330 
1.5880425E-03 
-0.1335113 


PHASE SHIFT ANALYSIS 


FREQ = 
0.8800000 
0.8800000 


w 
kp= lfreq 
AMPLITUDE; clamp, cmamp : 
joscil = Oitrans 
PHASE; сір, стр : 
AVERAGE DRAG, TOTAL DRAG : 
ETAS, WBAR : 


^ 


FREO SWEEP 
FREQ = 0.890000 


PHASE SHIF 
FREQ = 
0.8900000 
0.8900000 


ы 

kp= {freq 

AMPLITUDE; clamp, cmamp : 
foscil = Oitrans 
PHASE; сір; стр $ 
AVERAGE DRAG, TOTAL DRAG : 
ETAS, WBAR ° 


PHASE SHIF 

Š FREO 
0.8900000 

0.8900000 


w 
kp- ifreq 


AMPLITUDE; clamp, cmamp : 


Figure 3.2c UPOTFLUT 


0.8800000 
4 


0.2856959 

- 1 
271.3546 

-6.5488103Е-03 
0.1042214 


T ANALYSIS 
0.8900000 


5 
0.2319494 
= 0 
204.3115 


1.5859993Е-03 
=0.1304959 


T ANALYSIS 
0.8900000 
5 


0.2883391 


example 


2.9765518E-02 
-64.02047 


0.1048966 
-1.1639614E-02 


2.2874046Е-02 
36.47753 


-0.4246101 
-0.1231114 


3.0037180Е-02 
-64.15525 


0.1048108 
-1.1894441Е-02 


2. 325 В233Е-02 
36.17284 


-0.4322215 
-0.1256711 


3.0311935Е-02 
-64.29196 


0.1046759 
-1.2153637E-02 


2.3644408E-02 





output file 
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iloscil - Oitrans = 

PHASE; сір, стр : 271.6436 35.87401 
AVERAGE DRAG, TOTAL DRAG : -6.6649993Е-03 -0.4398900 
ETAS, WBAR : 0.1039310 -0.1282581 


FREQ SWEEP 
FREQ = 0.900000 


PHASE SHIFT ANALYSIS 
FREQ = 0.9000000 
w 0.9000000 
kp= 0.9000000 i freq 6 


AMPLITUDE; clamp, cmamp : 0.7320280 3.0592278Е-02 
1903611 = Oitrans = 0 

PHASE; clp, cmp : 204.7803 -64.43259 
AVERAGE DRAG, TOTAL DRAG : 1.5828811Е-03 0.1044701 
ETAS, WBAR : -0,1274617 -1.2418482Е-0?2 


PHASE SHIFT ANALYSIS 
FREQ = 0.9000000 
м 0.9000000 
kp= 0.9000000 1 freq 6 


AMPLITUDE; clamp, cmamp : 0.2909827 2.4039967E-02 
foscil = Oitrans = 1 

PHASE; clp, cmp 3 271.9288 35.58691 
AVERAGE DRAG, TOTAL DRAG : -6.7821071Е-03 -0.4476191 
ETAS, WBAR : 0.1036481 -0.1308680 


FREQ SWEEP 
FREQ - 0.910000 


PHASE SHIFT ANALYSIS 
FREQ - 0.9100000 
ы 0.9100000 
0.9100000 i freq 7 


AMPLITUDE; clamp, cmamp : 
foscil = Oitrans 


0.2321466 2:0865142Б-02 
= 0 

PHASE; € Lp; cmp 205,2529 -64.6122) 
AVERAGE DRAG, TOTAL DRAG 1.5800473Е-03 0.1042831 
ETAS, WBAR 2 -0-1245561 -1.2687406Е-02 


^ 


FHASE SHIFT ANALYSIS 
FREQ - 0.9100000 
ы 0.9100000 
Кр- 0.9100000 ifreq 7 


AMPLITUDE; clamp, cmamp : 0.2936268 2.4440434E-02 
foscil = Oitrans = 1 

PHASE; clp, cmp 212.2159 35.29980 
AVERAGE DRAG, TOTAL DRAG -6.9000078E-03 -0.4554005 
ETAS, WBAR 2 0.1033707 -0.1335003 


FREQ SWEEP 
FREQ = 0.920000 


PHASE SHIFT ANALYSIS 
ó FREQ = 0.9200000 

w 0.9200000 
Кр- 0.9200000 i freq 8 





Figure 3.2d UPOTFLUT example output file 
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A se wv eis 


AMPLITUDE: clamp, статр : 0.2 122 159 1.1118692ғҒ-02 


loscll - Oltrans - 0 

PHASE; clp, стр : 205.7373 -64.69626 
AVERACE DRAG, TOTAL DRAG : 1.5782372F-02 0.1041637 
FTAS, WRAR z -0.17187623 - 1. 2949491۴-07 


FHASE SHIFT ANALYSIS 


FRFO - 0.9200000 
ы 0.9200000 
kp- 0.9200000 (freet в 
AMPLITUDE: clamp, cmamp : 0.2962705 2.4841581Ғ-02 
loscil = Oltrans - 1 
PHASE; сір, стр : 272.5011 35.01268 
AVERAGE DRAC, TOTAL DRAG : -7.0186030Ғ-07 -0.4632331 
FIAS, WRAR : 0.101098? -0. 1361546 
FREQ SwrFP 
rRFO - 0.930000 


FHASE SHIFT ANALYSIS 


FREQ >= 0.9299999 
ы 0.9299999 
kp- 0.9299999 {freq 9 
AMPLITUDE; сізятр, статр : 0.212454Я 3.1417110Е-02 
loscil = Oltrans - 0 
РПА5Е; сір, стр : 206.2080 -61,02113 
AVERAGE DRAG, TOTAL DRAG : 1.5755624F-93 0.1039071 
ЕТЛ5, WRAR +: -0.1191657 -1.12216076-0? 
FHASE SWIFT ANALYSIS 
TRFO = 0.9799999 
w 0.9299999 
kp- 0.9799999 lfreg 9 
AMPLITUOF; clamp, cramp 1 0.298917] 2.8250112Ғ-0?2 
foscl) - Oltrans - 1 
PASE; clp, стр : 272.190] 31.172948 
AVERAGE DRAG, TOTAL DRAC : -7.13021086-03 -0.4711219 
ETAS, МАЛА : 0.1028326 -0.1 3882117 


FREQ SWEFP 
FREQ = 0.940000 


FHASE SHIFT ANALYSIS 


РАРО - 0,9199999 
w 0.9399999 
kr: 0.9199999 1frea 10 
ANELITUOF; clamp, cramp : 0.2326474 1.1699076ғ.-02 
1овс11 + Oltrans = 0 
PHASE; сір, cmp : 206.6846 -64.940821 
AVERAGE DRAG, TOTAL DRAG : 1.5715322f-03 0,1037211 
FTAS, WAAR 2 -0,1164354 -1.3496989Е-07 
а PHASE SHIFT ANALYSIS 
FRO - 0.9199999 
м 0.9199999 
kp: 0.9199999 {freq 10 
AMELITUDE: clamp, cmamp : ОП. 1015640 2.5659967ғ-02 
Тозс 11 = Oltrans > 1 
PHASE; сір, стр : 273.0792 31.4482) 
AVERACE ORAG, TOVAL DRAC : -1.2587]137F-03 -0.4790751 
E ҒІЛ5, МРАК : 0.1025746 -0.1415 104 
Number of Kp Values = 10 
Plunge Value/full Cherd (h/?bl - 0.0833 
Alpha ~ 3.0000 
PIVOT FOINT (a) Of Flastic Axis o -0.4000 
Half Chord (b) . 6.0000 
Alpha Uncoupled Nat. Freq. e %0.0000 
Plunge Uncoupled Nat. Freq . 72.0000 
1 alphe z 4.8410 
S$ alpha ы 0.6455 
Маза з 0.5379 
Alr Density ~ 0.0024 
Mass Ratio s 2.0000 
Dist. To Wing CG Aft of Elas. ^x1ss 0.20^0 


es EES ss KS: ern Ree e ës mr e mm mme re e E e 


Kp crit =- 0.9000000 


DIFF = 1.1507273Е-03 
SORTX = 0.9579428 
U crit > 1252. 694 


Figure 3.2e UPOTFLUT example output file 
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j. PLUNGE.in This file contains K,, Cy, RE, Cy, IM, Cmn RE, 
Суы ІМ 


k. FLUTPLOT.d This file contains K,, SORT (x) Re, SQRT (X) RE, 
SORT(X) IM 


3. VALIDATION 
The program was tested against some sample cases to 
check for code validity. The first case was taken from 
reference 6 example H1, p. 236. Figures 3.3 and 3.4 show 
plots of the FLUTPOT.d file. Figure 3.3 shows the initial 
look over a wide range of K, and after finding the approximate 
flutter location Figure 3.4 shows a closer look at the Kj, 


range of interest. This example calculated a U 1 of 


critica 
161.985 ft/sec. which compares favorably to the example value 
of 162 ft/sec. The next example was taken from NACA TR-685 
[ref.8] case #1 p. 8. Figures 3.5 and 3.6 again show the 
initial and final looks for this analysis. The example called 


for a U , ОЕ 567 miles/hr and the program returned a 


critica 
value of 570 miles/hr. Next, the code was tested over a range 
of w/w, ratios as done іп NACA ТЕ-685, р11., graph I-A(a). 


Figure 3.7 shows the comparison between the two methods. 


101 


090 СС () 





с90' 06" овас“ 09242%, 


6 


Q3DC-d5 
OCH 


1000МОУК) 119 


e 


СРО ОР О 


295 әЈашехя аоо[п]4 


(X)LYOS 


example 1 


Figure 3.3 Initial look at flutter, 
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Figure 3.4 Final look at flutter, 
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Figure 3.5 Initial look at flutter, example 2 
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Figure 3.6 Final look at flutter, 
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Mechanism of Flutter . The affects of changing tha ratio of plunge and pitch natural frequencies. 
Panal coda vse Theo. | (NACAOOO07, БОрапећ, Здаџ, ћ/25 =.0833, Зсусббсајс., ) 
| |(Mass= .53789 slugs/ft, lalpha = 4.84102, Salpha = .645468,b = ft. std day air dansity) 


Wh/W alpha Panel Ucrit.(ft/sec) Panel (-Ucnt/bWalpha) Theo (-Ucrit/bW alpha) 96Diff wrt Theo. 


| |. OO111| 896.6013! -1.66 -1.7 
ЕЕ. 1102 926.5568 -1.713994 -1.72' 


2.35% 
0.35 % 


[04 _ 9603143 -1.7783598 1.8) 1.20% 
| 405 994.4232 -1.84152 -1.9 3.08% 


[06 1036.167 -2.04 
0.7 : 





Mechanism of Flutter 
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Figure 3.7 w,/w, Calculations 


106 


5.94 96 





IV. FLOW VISUALIZATION EXPERIMENT 


A. INTRODUCTION 

The purpose of this experiment was to document the 
pasduction of thrust by a plunging airfoil. This was a 
preliminary experiment to better understand the vortex pattern 
produced by a plunging airfoil, and to examine the production 
of thrust using smoke flow visualization techniques. 

An explanation of what constitutes a propulsive vortical 
signature along with smoke flow visualization of the 
propulsive vortical patterns is given in Reference 7. In this 
reference, the explanation is given by contrasting the 
vortical pattern produced by a cylinder (drag) with the 
vortical pattern produced by a plunging airfoil (thrust). The 
cylinder produced a vortical sheet where the top row of 
vortices rotated clockwise and the bottom row of vortices 
rotated counterclockwise. This pattern induces a velocity 
component in the upstream direction (Biot-Savart law). In 
contrast, the plunging airfoil produced a clockwise rotating 
vortex sheet on the bottom row. This pattern induces a 
velocity component in the downstream direction. Reproduction 
of the flow visualization data from Reference 7 is shown in 


Figure 4.1. 
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Figure 4.1 View of flow over cylinder (top) and plunging 
airfoil (bottom) [Ref. 7] 
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В. THEORY 

A comparison was done using the incompressible panel code, 
U2DIIF. The purpose of this study was to examine the vortical 
pattern produced by the panel code, and determine if the 
vortical signature matched experimental results. The input to 
the panel code was set up to best match the conditions of the 
experiment described in the next section. The panel code was 
run using a plunge amplitude, h/2b equal to .1977, a reduced 
frequency of 1.8 and a zero mean AOA. The results of the 
vortical pattern are shown in Figure 4.2. Aside from the 
starting vortex, this is clearly a thrust producing vortical 
Sheet. Furthermore, the vortical pattern is similar to that 


produced by the experiment shown in Figures 4.10 and 4.11. 


C. EXPERIMENTAL SETUP 
1. Plunging Airfoil 

The plunging airfoil used in this paper was originally 
a wing taken from the rotor of a model helicopter. The wing 
was attached to a MB250 Shaker Table as shown in Figure 4.3 
The wing was made from a NACAO007 airfoil section and 
consisted of a 2.45" chord and a 22" span. The wing was built 
from a foam core and finished with a layer of graphite epoxy 
composite for added fatigue strength. The airfoil’s drive 


mechanism was a MB 250 Shaker Table capable of 1" total 
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Figure 4.2 Wake pattern produced by U2DIIF code 


110 


х 
мы 


хо SNS 
Š У 


е 
Ze 
san 


ЖЖ 


Ж 





Figure 4.3 Shaker table with wing 
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deflection. The Shaker Table was limited only by resonance 
frequencies of the wing which occurred around 20 Hz or 1200 
rpm. 
2. WIND TUNNEL 

The wind tunnel used in this experiment was a very low 
speed, low turbulence smoke tunnel. It is made of plexiglass 
walls and a contraction ratio of 2.8:1. The motor provides 
wind tunnel velocities between 0 and 10 feet per second (fps). 
The smoke was created using a Rosco smoke generator and piped 
into the tunnel in the test section using a small seven tube 
smoke rake constructed for this experiment. Figure 4.4 isa 
photograph of the wind tunnel and smoke rake used in this 


experiment. 


D. TEST PROCEDURE 

Testing was conducted in the low speed smoke tunnel under 
several different conditions. The speeds of the tunnel were 
approximately 1.04 fps, 1.47 fps, and 1.56 fps (measured 
visually). These low speeds allowed good pictures and the 
ability to get higher reduced frequencies without calling for 
too high a load on the wing. The actual plunging harmonic 
frequencies ranged from 1 to 15 Hz and amplitudes from 1/16" 
to 1" peak to peak. Тһе tunnel was initially turned off and 
the Shaker Table turned on with stagnant smoke in the tunnel. 
The purpose was to see if the 


plunging airfoil would draw the smoke through the tunnel like 
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Figure 4.4 Tunnel and smoke rake 


Pro 


a fan, thus showing the production of thrust by the plunging 
airfoil. 

Photos were taken using a Nikon 35mm camera and Kodak 
TMAX-400 ASA black and white film. The shutter speed was set 
to 1/125 seconds with an aperture setting of 4.0 for the light 
conditions. Film developing time was optimized at 9 minutes at 


75 degrees F. 


E. RESULTS AND DISCUSSION 

The result for the tunnel off condition flow visualization 
experiment was as expected. The wing in fact accelerated the 
smoke in itS vicinity. 

The result of the additional rake flow visualization 
experiments are shown in Figures 4.5-4.14. Figure 4.5 shows 
the stationary airfoil at zero degree AOA. The Reynolds number 
(based on airfoil chord) is 10,000. It can be seen that the 
airfoil produces a small wake with the boundary layer mostly 
attached. Figures 4.6 through 4.14 show the vortical wake 
flow patterns produced by plunge oscillations at various 
frequencies as indicated. Most of these pictures reveal the 
propulsive vortical street pattern discovered in Reference 7. 
Previous experiments by Neace, [ref.9] found that the tunnel 
was too small for the airfoil size used, but the airfoil size 
for the present experiment seemed to be optimum, as seen by 
the long trail of vortices. The vortical patterns show that 


the bottom vortex is rotating clockwise, and the top vortex is 
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Figure 4.5 Steady airfoii 1.56 ft 
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Figure 4.6 Кр 1.8008, h/2b-.1977, 1.56 ft/s 


хе хок ко» 











2 
Е: 
4 
E 
| 








56 ft/s 


1 


h/2b=.10204, 
TI 


2467; 


- 2 


Figure 4.7 Кр 
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Figure 4.9 Ko = 4.112, h/2b=.14031, 1.56 ft/s 
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Figure 4.12 Кр = 7.401, h/2b=.03826, 1.56 ft/s 
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rotating counterclockwise, which is а thrust producing 
vortical sheet. It can be seen in the pictures that the 
frequency greatly affects the vortical strength (size). 


Increasing the frequency leads to an increase in wake 


BEE 
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У. LIFT ENHANCEMENT PRODUCED BY A PLUNGING AIRFOIL 


A. THEORY 

Chapter IV demonstrated the propulsive capability of a 
plunging airfoil. The production of thrust implies the 
generation of a jet flow which, in turn, may be utilized asa 
boundary layer control device. Therefore, an additional test 
was conducted in the NPS smoke tunnel in order to explore the 


feasibility of this concept. 


B. SETUP 

The same NACA0007 plunging wing was used as in chapter IV 
with a different driving mechanism. The wing was mounted to 
an ELECTRO-SEIS Model 113 Shaker Table by APS Dynamics, Inc. 
The shaker was located below the test section of the NPS smoke 
tunnel (Figure 5.1). The plunging airfoil was mounted to 
struts at both ends to prevent excessive bending while 
plunging. The large airfoil is a cambered profile taken from 
the rotor of a full size helicopter (13" chord, 2" Thickness, 
modified NACA airfoil, Reynolds Number of 52,000). The wing is 
Suspended from the tunnel ceiling as shown in the flow 
pictures. The design allowed for full movement of the big wing 


to position it in the vicinity of ‘the plunging alxromm 
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С. WIND TUNNEL 

This study used the Naval Postgraduate School’s flow 
visualization wind tunnel. The tunnel iS an open-circuit one, 
with air entering an inlet that measures 4.5 m X 4.5 m 
(15'X15'). As the air enters the tunnel, it passes through a 
7.5-cm long honeycomb. A 9:1 ratio Square contraction cone 
directs the flow into a test section that is 1.5 m X 1.5 m 
(5'X5'), and 6.7 m (22’) long, as seen in Figure 5.2. The 
flow is then exhausted into the atmosphere through a fan, 
which uses variable pitch blades to control the speed of the 
flow. The speed control toggle switch is located right below 
the red and green on/off switch located in the left side of 
the tunnel control room. The tunnel speed was determined using 
a digital manometer which was verified for accuracy (Figure 
E 

An E is located on the side of the tunnel. 
A glass window, 1.6 ш Х 1.1 ш (5.2” X 3.4”), provides the 
primary viewing area from the observation room and a second 
one, 0.4 m X 1.23 m (1.33'X4'), is located in the tunnel’s 
roof. The main viewing window had sufficient area for most of 
the photography, with the top window used for illumination. A 
circular turntable was located on the floor of the test 
section [ref.11] which allowed for easy access to the shaker 
table. The walls and floor of the test section were flat black 


for low light reflectivity. 
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Figure 5.1 Shaker table setup below tunnel 
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Figure 5.2 Tunnel Layout 
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Figure 5.3 Tunnel speed check 
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D. SMOKE GENERATION 

The smoke was generated in the Rosco smoke/fog machine. 
Many different smoke injection techniques were tried but with 
less than satisfactory results. Smoke rakes were first tried 
outside the tunnel with the tube number varying from 2 to 30 
tubes. The tubes were inserted in the honeycomb and also 
separated different distances from the inlet of the tunnel. 
The rake was also tried inside the test section with very bad 
results (smoke dispersed immediately). Problems ranged from 
lack of smoke and turbulence when enough smoke was present. 
The Rosco machine at its lowest setting was producing a very 
high smoke volume and whenever the flow was restricted by a 
smoke rake the smoke production went way down. The final 
technique adopted was very simple. The smoke output was sent 
directly from the machine to a 1" nozzle which was manually 
waved at the entrance of the tunnel to make a steady cloud. 
The steady cloud was gradually pulled into the tunnel, 


producing a thick smoke sheet in the test section. 


E. PHOTOGRAPHY 

Photos were taken using a Nikon 35mm camera and Kodak 
TMAX-400 ASA black and white film. Тһе film speed was set to 
1/250 seconds with an aperture setting of 4.0 for the light 
conditions. Film developing time was optimized at 9 minutes at 


75 degrees F. 


TT 


Е. EXPERIMENTAL PROCEDURES 

The first step involved a look at the large airfoil to 
verify normal flow patterns (Figure 5.4 and 5.5) and find the 
AOA for initial trailing edge separation (Figure 5.6). Next, 
the plunging airfoil was placed in the tunnel by itself anda 
run was made to verify the propulsive capability in the larger 
tunnel at higher speeds. As seen in Figure 5.7, the airfoil 
produced a drag vortical flow. Figure 5.8 and 5.9 shows the 
propulsive pattern of the propulsive airfoil. Finally, the two 
airfoils were placed in close relative position to see the 


interference effect between the two airfoils. 


G. RESULTS AND DISCUSSION 

Several airfoil position combinations were studied, as 
shown in Figures 5.10 through 5.15. Figure 5.10 and 5.11 show 
the plunging airfoil located at approximately .65 chord of the 
large airfoil. Figure 5.12 through 5.15 show the plunging 
airfoil located at approximately .75 chord. 

The differences between the plunging on and off condition 
were not easy to see with the eye but pictures indeed showed 
some differences between the two conditions. A shortcoming of 
this experiment was the inability of the plunging airfoil to 
run parallel with the large airfoil. Additionally, sizing of 
and relative positioning of the two airfoils was not optimized 
to give best results. The two airfoils were chosen from the 


resources available and time constraints prevented a more 
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Figure 5.4 Large airfoil at zero AOA 
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5 Large airfoil at 10 degrees AOA 


Figure 5. 
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Figure 5.6 Large airfoil at 12 degrees AOA (stall) 
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Figure 5.7 Plunge airfoil 








Figure 5.8 Plunging airfoil K =1.71 
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Figure 5.10 Large airfoil at 12 degrees AOA, steady plunging 
airfoil position 1 


хх N 


Figure 5.11 Large airfoil at 12 degrees 
airfoil К,-3.42, position 1 
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Figure 5.12 Large airfoil at 12 degrees AOA, plunging 
airfoil steady, position 2 


ка 
d 
+? 


1) 


BR i П 





plunging 


142 


3.42, 


Figure 5.13 Large airfoil at 12 degrees AOA, 
foil position 2 
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Figure 5.15 Large airfoil at 14 degrees AOA, plunging 
airfoil Кр=3.42, position 2 


144 


detailed investigation of the interference effects between the 


two airfoils. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. SINGLE AIRFOIL ANALYSIS 

The modified version of U2DIIF (UPOT) can perform 
aerodynamic calculations over any range of reduced 
frequencies. The nonlinear theory presented here for harmonic 
motion, and the phase relationships that exist between the 
airfoil motion and the aerodynamic forces have been 
extensively verified by comparison with Theodorsen’s linear 
theory. Furthermore, this panel code was applied to the 
analysis of incompressible bending-torsion airfoil flutter. 
Again, excellent agreement with the classical Theodorsen 
analysis was obtained. 

Access to faster computational means is recommended to 
абд е КС time needed to predict the flutter points. The 
code should be modified to incorporate three-dimensional 


calculations which would help solve more difficult flutter 


problems. 


B. FLOW VISUALIZATION EXPERIMENTS 

The flow visualization experiment successfully showed the 
development of thrust produced by a plunging airfoil. The 
enhanced lift experiment, on the other hand, was not a 
complete success. The smoke visualization presented 


difficulties that were not satisfactorily overcome. As a 
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result, the pictures taken were somewhat inconclusive. 
Futhermore, the angle of attack of the oscillating airfoil 
could not be changed thus making it difficult to achieve a 
flow condition conducive to lift enhancement. 

It is recommended that further experiments be conducted in 
the low speed smoke tunnel with a shaker table capable of 
moving an airfoil at harmonic frequencies near 40 HZ. 
Additionally, the airfoil must be modified to allow change of 
AOA. Finally, the Rosco smoke machine output volume must be 
modified to permit much lower smoke output. This final point 
proved to be the single largest detriment to the visualization 


experiment. 
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